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Abstract 

We present a measurement of the top quark mass using a sample of ti decays 
into an electron or a muon, a neutrino, and four jets. The data were collected 
in pp collisions at yfs = 1.8 TeV with the Collider Detector at Fermilab and 
correspond to an integrated luminosity of 109 pb _1 . We measure the top 
quark mass to be 175.9 ± 4.8(stat.) ± 4.9(syst.) GeV/c 2 . 

14.65.Ha, 13.85.Qk, 13.85.Ni 
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The top quark mass is a fundamental parameter of the standard model and is needed 
for extracting other parameters from precision electroweak measurements. The first direct 
measurement of its value was made by CDF |T[ and was based on 19 pb _1 of data. Updated 
measurements were reported by both the CDF and D0 collaborations using significantly 
more data In this paper we present a new measurement of the top quark mass with 

greatly improved precision, using our entire data sample from the 1992-1995 runs, which 
corresponds to a total integrated luminosity of 109 ± 7 pb -1 J7|. This new measurement 
supersedes the results reported in |l]fj]. 

Within the standard model, the top quark decays more than 99% of the time into Wb. 
The W boson can then decay to a quark- ant iquark or lepton-neutrino pair. The measurement 
presented here uses events with a tt pair decaying in the "lepton+jets" channel. This channel 
is characterized by a single high-P T |I[ lepton (electron or muon) and missing transverse 
energy from a W — > iv decay, plus several jets coming from a hadronically decaying W 
boson and from the b quarks from the top quark decays. Jets formed by the fragmentation 
of b quarks can be identified ("tagged") either by reconstructing secondary vertices from 
b hadron decays with the silicon vertex detector (SVX tagging), or by finding additional 
leptons from semileptonic b decays (SLT tagging). The SVX and SLT tagging algorithms 
are described in Ref. M. 

To be used for the mass measurement, events must contain a single isolated electron 
(muon) with E T (P T ) > 20 GeV (GeV/c) in the central region of the detector (\r)\ < 1) and 
missing transverse energy, t T > 20 GeV, indicating the presence of a neutrino. At least four 
jets are required in each event, three of which must have an observed E T > 15 GeV and 
\q\ < 2. In order to increase the acceptance, we relax the requirements on the fourth jet to 
be E T > 8 GeV and \r}\ < 2.4, provided one of the four leading jets is tagged by the SVX 
or SLT algorithms. SVX tags are only allowed on jets with observed E T > 15 GeV, while 
SLT tags are allowed on jets with E T > 8 GeV. If no such tag is present, the fourth jet must 
satisfy the same E T and rj requirements as the first three. All jets in this analysis are formed 
as clusters of calorimeter towers within cones of fixed radius AR = y/ Ai] 2 + A0 2 = 0.4 ||. 
The above selection defines our mass sample, which contains 83 events. 

Measurement of the top quark mass begins by fitting each event in the sample to the 
hypothesis of tt production followed by decay in the lepton+jets channel: 

pp — ► tt + X 

t >W + b >i + vb qq'b 

t >W~l >qq'b ° r i-vb. 

The 3-momenta of the lepton and the b, b, q and q' quarks are measured from the observed 
lepton and four leading jets in the event; the mass of the b is set to 5 GeV/c 2 , that of q 
and q 1 to 0.5 GeV/c 2 . The neutrino mass is assumed to be zero and its momentum is not 
measured, thereby yielding three unknowns. The two transverse momentum components of 
X are measured from the extra jets in the event and the energy that is detected but not 
collected in jet or electron clusters. Five constraints are applied: the transverse momentum 
components of the entire tt + X system must be zero, the invariant masses of the lepton- 
neutrino and q-q' pairs must each equal the W boson mass, and the mass of the top quark 
must equal that of the antitop quark. The problem therefore has two extra constraints 
and is solved by a standard ^-minimization technique. The output of each event fit is a 
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reconstructed top mass M rec and a \ 2 value quantifying how well the event is described by 
the tt hypothesis. 

Electron energies and muon momenta entering the fit are measured with the calorimeter 
and tracking chambers, respectively ||. Jet energies are corrected for losses in cracks be- 
tween detector components, absolute energy scale, contributions from the underlying event 
and multiple interactions, and losses outside the clustering cone. These corrections are de- 
termined from a combination of Monte Carlo simulations and data [Iu|. The four leading 



jets in a tt candidate event undergo an additional energy correction that depends on the 
type of parton they are assigned to in the fit: a light quark, a hadronically decaying b quark, 
or a b quark that decayed semileptonically [|TJ. This parton-specific correction was derived 
from a study of tt events generated with the herwig Monte Carlo program [|Tl| , [T^1 . 

There are twelve distinct ways of assigning the four leading jets to the four partons 
b, b, q, and q' . In addition, there is a quadratic ambiguity in the determination of the 
longitudinal component of the neutrino momentum. This yields up to twenty-four different 
configurations for reconstructing an event according to the tt hypothesis. We require that 
SVX or SLT-tagged jets be assigned to 6-partons and choose the configuration with lowest 
X 2 - Events with \ 2 > 10 are rejected. In the mass sample, 76 out of 83 events remain 
after this cut. When all parton-jet assignments are correctly made, the resolution of the 
reconstructed mass is 13 GeV/c 2 for a top mass of 175 GeV/c 2 . 

A maximum-likelihood method is used to extract a top mass measurement from a sample 
of events which have been reconstructed according to the tt hypothesis. An essential ingredi- 
ent of the likelihood function is the probability density / s (M rec |M top ) to reconstruct a mass 
M rec from a tt event if the true top mass is M top . In past publications [lip] we estimated f s 
for a discrete set of M top values by smoothing histograms of M rec for events from a herwig 
Monte Carlo calculation. In the present analysis we parameterize f s as a smooth function 
of both M rec and M top [0. This new approach yields a consistent, M top -dependent way 
of dealing with low statistics in the tails of the M rec histograms and produces a continuous 
likelihood shape from which the top mass and its uncertainty can be extracted. The proba- 
bility density /&(M rec ) for reconstructing a mass M rec from a background event is obtained 



by fitting a smooth function to a mass distribution generated with the VECBOS |b|] W+jets 
Monte Carlo program. 

The likelihood function is the product of three factors: 

£ £ shape X Csbackgr X Csparam j (1) 

where C s ha P e represents the joint probability density for a sample of N reconstructed masses 
Mi to be drawn from a population with a background fraction x&: 

N 

Cshape = II K 1 - X b) fs{Mi\M top ) + X b f b {Mi)} . 
i = l 

The fraction Xb is constrained by an independent measurement that is summarized by the 
background likelihood Cbackgr- The function C param allows the parameterizations of f s and 
to vary within the uncertainties returned by the fits to the herwig and VECBOS histograms 
of M rec . By including C param in the likelihood definition, the uncertainty due to the finite 
statistics of these histograms is incorporated into the statistical uncertainty on the measured 
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top mass. The likelihood C is maximized with respect to M top , x&, and the parameters that 
define the shapes of f s and /&. 

The precision of the top quark mass measurement is expected to increase with the 
number of observed events, the signal-over-background ratio, and the narrowness of the 
reconstructed-mass distribution. These characteristics vary significantly between samples 
with different b tagging requirements. Therefore, to make optimal use of all the available in- 
formation, we partition the mass sample into non-overlapping subsamples, define subsample 
likelihoods according to eq. (p]), and maximize the product of these likelihoods to determine 
the top mass and its uncertainty JT5| . The use of non-overlapping subsamples ensures that 
the corresponding likelihoods are statistically uncorrelated. Monte Carlo studies show that 
an optimum partition is made up of four subsamples: events with a single SVX tag, events 
with two SVX tags, events with an SLT tag but no SVX tag, and events with no tag but 
with the tighter kinematic requirement of four jets with E T > 15 GeV and \t)\ < 2. 

The calculation of the expected background content of each subsample starts from the 
background calculation performed on the W+ > 3-jet sample for the ti cross section mea- 
surement 0. The extrapolation to the mass subsamples takes into account the additional 
requirement of a fourth jet, the \ 2 cu ^ on event reconstruction, and the fact that SVX and 
SLT tags are only counted if they are on one of the four leading jets. The efficiencies of 
these requirements are determined from Monte Carlo studies. They are used together with 
background rates and tagging efficiencies from the cross section analysis to predict the to- 
tal number of events in each mass subsample as a function of the unknown numbers of ti 
and VF+jet events in the combined sample. These unknowns are estimated by maximiz- 
ing a multinomial likelihood that constrains the predicted subsample sizes to the observed 
ones. This procedure generates the expected background fractions shown in Table | and the 
background likelihood Cbackgr used in eq. (JTJ). 

Approximately 67% of the background in the entire mass sample comes from W+jet 
events. Another 20% consists of multijet events where a jet is misidentified as a lepton and 
bb events with a b hadron decaying semileptonically. The remaining 13% is made up of Z+jet 
events where the Z-boson decays leptonically, events with a WW, WZ or ZZ diboson, and 
single-top production. We have compared the reconstructed-mass distributions in vecbos 
and data for three event selections that are expected to be depleted in ti events [ 16] . These 
selections are slight variations of the mass sample selection. The first one requires that the 
primary lepton be an electron with a pseudo-rapidity in the range 1.1 < \r)\ < 2.4 instead 
of \rj\ < 1, and yields 26 data events. The second one requires at least four jets with E T > 8 
GeV and \rj\ < 2.4, but no more than two jets with E T > 15 GeV and \rj\ < 2. This results 
in 243 data events. The third selection requires events with a non-isolated primary lepton 
and yields 164 data events. In all three Kolmogorov-Smirnov test applied to the 

comparison of VECBOS and data yields a confidence level of at least 30%. We therefore use 
the vecbos calculation to determine the shape of fb for the likelihood function. 

The reconstructed-mass distribution of the sum of the four subsamples is plotted in 
Figure |I|. The inset shows the shape of the corresponding sum of negative log-likelihoods 
as a function of top mass. From this we measure M top = 175.9 ± 4.8 GeV/c 2 , where the 
uncertainty corresponds to a half-unit change in the negative log-likelihood with respect to 
its minimum. Monte Carlo studies on mass samples similar to ours yield an 11% probability 
for obtaining a statistical uncertainty of this size or smaller. The background fractions 
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returned by the fit agree with the x° numbers listed in Table |. To judge the goodness of 
the fit of the combined M rec distribution, we performed a Kolmogorov-Smirnov test and 
obtained a confidence level of 64%. The reconstructed-mass distribution in each of the four 
subsamples is compared to the result of the combined fit in Figure |2|. The insets show the 
results of likelihood fits performed separately in each of the four subsamples. The mass 
measurements obtained from these fits are consistent with each other, as shown in Table |. 

We list the systematic uncertainties in Table [IT]. The largest one comes from the jet 
energy measurement. Each of the jet energy corrections described earlier carries with it a 
separate, energy-dependent uncertainty (TD|. Recent studies of soft gluon radiation outside 



the jet clustering cone have reduced the uncertainty from this source to 2.5% for a jet 
with observed E T > 40 GeV. For an observed jet E T of 40 GeV, the total uncertainty on 
the corrected E T varies between 3.4 and 5.6% depending on the proximity of the jet to 
cracks between detector components. We have checked the jet correction procedure and the 
evaluation of the jet energy scale uncertainty with events containing a leptonically decaying 
Z boson and one jet. A study of how the transverse momentum of the jet balances that 
of the Z decay products finds that the observed ratio of [P T (Z) — P T (jet)]/P T (Z) differs 
by 3.2 ± 1.5(stat.) ± 4.1(syst.) % from Monte Carlo simulations. The 4.1% systematic 
uncertainty is due to the jet energy scale only. Since the difference is consistent with zero, 
this study independently confirms the soundness of our estimate of the jet energy scale 
uncertainty. A further confirmation was obtained by measuring the mass of the W boson 
from its hadronic decay modes, using a sample of ti candidate events in the lepton+jets 
channel. This measurement yields 77.2 ± 3.5(stat.) ± 2.9(syst.) GeV/c 2 |17 . 



The second largest systematic uncertainty is due to high transverse momentum gluons 
that are radiated from the initial or final state of a ti event and sometimes take the place of 
a ti decay product among the four leading jets. This uncertainty was determined with the 
pythia Monte Carlo calculation Jl8| by separately studying the effect of extra jets coming 
from initial and final state radiation. 

The uncertainty in the modeling of the background mass distribution was estimated by 
varying the Q 2 scale in vecbos. Additional sources of uncertainty include the kinematical 



bias introduced by b tagging and the choice of parton distribution functions (cteq4l (19 
vs. mrsdO'). The sum in quadrature of all the systematic uncertainties is 4.9 GeV/c 2 . We 
have investigated the effect of using Monte Carlo calculations other than herwig to model 
ti events. Whereas pythia yields the same measured mass, isajet [[HJ leads to a +1.5 
GeV/c 2 shift. We do not include this as a separate uncertainty since the main difference 
between these calculations, namely the modeling of gluon radiation and jet fragmentation, 
is already accounted for in our analysis of other systematic uncertainties. 

In summary, we have measured the top quark mass to be 175.9 ± 4.8(stat.) ± 4.9(syst.) 
GeV/c 2 . This is the most precise determination of the top mass to date. A new technique 
for optimizing the use of the information provided by the tagging algorithms has resulted 
in a smaller statistical uncertainty, and a better understanding of the jet energy scale has 
led to a reduced systematic uncertainty. In addition, the probability densities for recon- 
structed masses are now fully parameterized, which simplifies the likelihood analysis and 
the treatment of the finite statistics of the Monte Carlo event samples. 

We thank the Fermilab staff and the technical staffs of the participating institutions for 
their vital contributions. This work is supported by the U.S. Department of Energy and 
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TABLES 



TABLE I. Subsamples of W+ > 4-jet events that are used for the top quark mass measurement. 
For each subsample, the number of observed events N f, s , the expected background fraction xP, and 
the measured top mass M top are shown. Uncertainties on the measured top mass are statistical 
only. 



Subsample 


N obs 


4 
(%) 


Measured M top 
(GeV/c 2 ) 


SVX double tag 
SVX single tag 
SLT tag (no SVX) 
No tag (E T (j 4 ) > 15 GeV) 


5 

15 
14 

42 


5±3 
13 ±5 
40 ±9 
56 ± 15 


170.1 ±9.3 
178.0 ±7.9 

142lf| 
181.0 ±9.0 


TABLE II. List of systematic uncertainties 


on the final top quark mass measurement. 


Source 






Value (GeV/c 2 ) 


Jet energy measurement 
Initial and final state radiation 
Shape of background spectrum 

6 tag bias 
Parton distribution functions 






4.4 
1.8 
1.3 
0.4 
0.3 


Total 






4.9 
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FIGURES 




FIG. 1. Reconstructed-mass distribution of the four mass subsamples combined. The data 
(points) are compared with the result of the combined fit (dark shading) and with the background 
component of the fit (light shading). The inset shows the variation of the combined negative 
log- likelihood with M top . 
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FIG. 2. Reconstructed-mass distributions in each of the four mass subsamples. Each plot 
shows the data (points), the result of the combined fit to top+background (dark shading), and the 
background component of the fit (light shading). The insets show the variation of the negative 
log-likelihoods with Mt op for the separate subsample fits. 
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